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Abstract 

We describe arithmetic computations in terms of operations on 
some well known free algebras (SIS, S2S and ordered rooted 
binary trees) while emphasizing the common structure present in 
all them when seen as isomorphic with the set of natural numbers. 

Constructors and deconstructors seen through an initial algebra 
semantics are generalized to recursively defined functions obeying 
similar laws. Implementations using Scala's apply and unapply 
are discussed together with an application to a realistic arbitrary 
size arithmetic package written in Scala, based on the free algebra 
of rooted ordered binary trees, which also supports rational number 
operations through an extension to signed rationals of the Calkin- 
Wilf bijection. 
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1. Introduction 

Classical mathematics frequently uses functions defined on equiv- 
alence classes (e.g. modular arithmetic, factor objects in algebraic 
structures) provided that it can prove that the choice of a represen- 
tative in the class is irrelevant. 

On the other hand, when working with proof assistants like 
Coq [9], based on type theory and its computationally refined 
extensions like the Calculus of Construction [4], one cannot avoid 
noticing the prevalence of data types corresponding to free objects, 
on top of which everything else is built in the form of canonical 
representations. 

Category-theory based descriptions of Peano arithmetic fit nat- 
urally in the general view that data types are initial algebras - in 
this case the initial algebra generated by the successor function, as 
a provider of the canonical representation of natural numbers. Of 
course, a critical element in choosing such free algebras is compu- 
tational efficiency of the operations one wants to perform on them, 
in terms of low time and space complexity. For instance Coq for- 
malizations of natural numbers typically use binary representations 
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while keeping the Peano arithmetic view when more convenient in 
proofs [9]. Note also that free algebras corresponding to one and 
two successor arithmetic (SIS and S2S) have been used as a basis 
for decidable weak arithmetic systems like [2] and [1 1]. It has been 
shown recently in [16, 17] that the initial algebra of ordered rooted 
binary trees corresponding to the language of Godel's System T 
types [7] can be used as a the language of arithmetic representa- 
tions, with hyper-exponential gains when handling numbers built 

2 

from "towers of exponents" like 2 2 . Independently, this view is 
confirmed by the suggestion to use A-terms as a form of universal 
data compression tool [8] and by deriving bijective encodings of 
data types using a game-based mechanism [18]. 

These results suggest a free algebra based reconstruction of 
fundamental data types that are relevant as building blocks for 
finite mathematics and computer science. We will sketch in this 
paper an (elementary, not involving category theory) foundation for 
arithmetic computations with free algebras, in which construction 
of sets, sequences, graphs, etc. can be further carried out along the 
lines of [14, 15, 17]. 

The paper is organized as follows. We define in section 2 iso- 
morphisms between the free algebras of signatures consisting of 
one constant and respectively, of one successor (S), two successors 
(0 and I) and a free magma constructor (C). 

To enable computations with the objects of the free algebras, 
we discuss the use of generalized constructors / destructors derived 
from these free algebras using the apply / unapply constructs avail- 
able in Scala (section 3). As an application, a complete arbitrary 
size rational arithmetic package using the Calkin- Wilf bijection be- 
tween positive rationals and natural numbers is described in section 
4. Sections 5 and 6 discuss related work and our conclusions. 

2. Free Algebras and Data Types 

DEFINITION 1 . Let a be a signature consisting of an alphabet of 
constants ( called generators) and an alphabet of function symbols 
(called constructors) with various arities. We define the free algebra 
A a of signature a inductively as the smallest set such that: 

1. if c is a constant of a then c £ Ac 

2. if f is an n-argument function symbol of a, then Vi, < i < 
U,ti E A a =>■ /(t , . . . ,U, . . .,t n -l) G A a . 

We will write c/0 for constants and f /n for function symbols with 
n arguments belonging to a given a signature. 

More general definitions, e.g. as initial objects in the category of 
algebraic structures, are also used in the literature and a close rela- 
tion exists with term algebras distinguishing between function con- 
structors (generating the Herbrand Universe) and predicate con- 
structors (generating the Herbrand Base). 

Recursive data types in programming languages like Haskell, 
ML, Scala can be seen as a notation for free algebras. We refer to 
[19] for a clear and convincing description of this connection. 

For instance, the Haskell declarations 



data AlgU = U I S AlgU 

data AlgB = B I AlgB I I AlgB 

data AlgT = T I C AlgT AlgT 

correspond, respectively to 

• the free algebra AlgU with a single generator U and unary 
constructor S (that can be seen as part of the language of Peano 
or Robinson arithmetic, or the decidable (W)S1S system, [2]) 

• the free algebra AlgB with single generator B and two unary 
constructors and I (corresponding to the language of the 
decidable system (W)S2S [11]), as well as "bijective base-2" 
number notation [20] 

• the free algebra AlgT with single generator T and one binary 
constructor C (essentially the same thing as the free magma 
generated by T). 

The set-theoretical construction corresponding to the " I " operation 
is disjoint union and the data types correspond to infinite sets gen- 
erated by applying the respective constructors repeatedly. The set- 
theoretical interpretation of "self-reference" in such data type defi- 
nitions can be seen as fixpoint operation on sets of natural numbers 
as shown in the Poj construction used by Dana Scott in defining the 
denotational semantics for various A-calculus constructs [13]. 

We will next "instantiate" some general results to make the un- 
derlying mathematics as elementary and self-contained as possible. 
While category theory is frequently used as the mathematical back- 
ing for data-types, we will provide here a simple set theory-based 
formalism, along the lines of [1]. 

We will start with the elementary mathematics behind the AlgT 
data type and follow with an outline for a similar treatment of AlgU 
and AlgB. 

2.1 The free magma of ordered rooted binary trees with 
empty leaves 

DEFINITION 2. A set M with a (total) binary operation * is called 
a magma. 

DEFINITION 3. A morphism between two magmas M and M' is a 
function f : M — >• M' such that f(x * y) = f(x) * f(y). 

Let X be a set. We define the sets M n (X) inductively as 
follows: M\(X) — X and for n > 1, M„(X) is the disjoint union 
of the sets M k {X) x M n - k (X) for < k < n. Let M(X) be 
the disjoint union of the family of sets M n (X) for n > 0. We 
identify each set M n with its canonical image in M(X). Then for 
w £ M n (X), we call n the length of w and denote it l(w). Let 
w, w' € M(X) and let p = l(w) and q — l(w'). The image of 
(w, w') £ M p x M q under the canonical injection in M(X) is 
called the composition of w and w' and is denoted w * w'. 

When X = {T} where T is interpreted as the "empty" leaf 
of ordered rooted binary trees, the elements of M n can be seen as 
ordered rooted binary trees with n leaves while the composition 
operation "*" represents joining two trees at their roots to form a 
new tree. 

DEFINITION 4. The set M(X) with the composition operation 
(w, w') — > w * w' is called the free magma generated by X. 

PROPOSITION 1. Let M be a magma. Then every mapping u : 
X — ► M can be extended in a unique way to a morphism ofM(X) 
into M. 

Proof We define inductively the mappings /„ : M m (X) — > M 
as follows: For n = l,/i = /. For n > l,Vp € {1, ..,n — 
1}, f n (w * w') = f p (w) * f„- p {w'). Let g : M(X) -> M such 
thatVn > 0,Va; € M n (X),g(x) = f n (x). Then g is the unique 
morphism of M(X) into M which extends f. 



Note that this property corresponds of the initial algebra [19] 
view of the corresponding (ordered, rooted) binary tree data type. 

Definition 5. If u : X -> Y, we denote M(u) : M(X) -> 
M(Y) the unique morphism of magmas defined by the construction 
in Proposition 1. 

\iv:Y—¥Z then the morphism M(v) o M(u) extends v o u : 
X -> Z and therefore M(v) o M(u) = M(v o u). 

PROPOSITION 2. If u : X — > Y is respectively injective, surjec- 
tive, bijective then so is M(u). 

It follows that 

Proposition 3. If X = {x} and Y = {y} and u : X -»• Y is 
the bijection such that f(x) = y, then M(u) : M(X) -> M(Y) 
is a bijective morphism (i.e. an isomorphism) of free magmas. 

Proof If X is empty so is M(X), hence u is injective. If u is in- 
jective, then 3u' : Y — > X, u o u = idu(x) where idu(x) 
denotes the identity mapping of M(X). Then M(u') o M(u) = 
M(u' o u) = idM(x) an d hence M(u) is injective. If u is surjec- 
tive, then 3-u' : 7 -> X,uou' = id M (y) ■ Then M(u) o M(u) = 
M(u o u') — id M (Y) and hence M(u) is surjective. If u is bijec- 
tive, than it is injective and surjective and so is M(u). 

We will identify the data type AlgT with the free magma gener- 
ated by the set {T} and denote its binary operation x * y as C x y. 
It corresponds to the free algebra (that we will also denote AlgT) 
defined by the signature {T/0 , C/2}. 

We can now instantiate the results described by the previous 
propositions to AlgT: 

PROPOSITION 4. Let X be an algebra defined by a constant t and 
a binary operation c. Then there's a unique morphism f : AlgT — > 
X that verifies 

f(T) = t (1) 

f(C(x,y))=c(f(x),f(y)) (2) 
Moreover, if X is a free algebra then f is an isomorphism. 

Proof It follows from Proposition 2 and equation /(T) = t, given 
that / is a bijection between the singleton sets {T} and {t}. 

2.2 The One Successor and Two Successors Free Algebras 

The one successor free algebra (also known as unary natural num- 
bers or Peano algebra, as well as the language of the monoid {0}* 
and the decidable systems WS1S and SIS) is defined by the sig- 
nature {U/0, S/l}, where U is a constant (seen as zero) and S is 
the unary successor function. We will denote AlgU this algebra and 
identify it with its corresponding data type. 

We state an analogue of Proposition 4 for the free algebra AlgU. 

PROPOSITION 5 . Let X be an algebra defined by a constant u and 
a unary operation s. Then there's a unique morphism f : AlgU — > 
X that verifies 

f(U) = u (3) 

f(S(x)) = s(f(x)) (4) 
Moreover, if X is a free algebra then f is an isomorphism. 

Note that following the usual identification of data types and initial 
algebras, AlgU corresponds to the initial algebra "1 + _" through 
the operation g =<U,S> seen as a bijection g : 1 + N — > N. 

The two successor free algebra (also known as bijective base- 
2 natural numbers or Peano algebra, as well as the language of 
the monoid {0, 1}* and the decidable systems WS2S and S2S ) is 



defined by the signature {B/0 , 0/1, 1/1} where B is a constant 
(seen as the empty sequence) and 0, I are two unary successor 
functions. We will denote AlgB this algebra and identify it with its 
corresponding data type. 

We can state an analogue of Proposition 4 for the free algebra 
AlgB. 

PROPOSITION 6. Let X be an algebra defined by a constant b 
and a two unary operations o, i. Then there 's a unique morphism 
f : AlgB — > X that verifies 

f(B) = b (5) 
f(0(x)) = o(f(x)) (6) 

/(/(i)) = i(f(x)) (7) 
Moreover, if X is a free algebra then f is an isomorphism. 

These observations suggest that for defining isomorphisms be- 
tween AlgU, AlgB and AlgT that enable a complete set of equiva- 
lent arithmetic (and later set-theoretic) operations on each of them, 
we will need a mechanism to prove such equivalences. To this end, 
it will be enough to prove that such non-constructor operations also 
form free algebras of matching signatures. 

We will call terms the elements of our initial algebras. 

3. Generalized Constructors 

The iso-functors supporting the equivalence between actual con- 
structors and their recursively defined function counterparts sug- 
gest exploring programming language constructs that treat them in 
a similar way. For instance it makes sense to extend "constructor- 
only benefits" like pattern matching to their function counterparts. 

Fortunately, constructors/deconstructors generalized to arbi- 
trary functions are available in Scala through apply /unapply 
methods and in Haskell through a special notation implement- 
ing views, under the implicit assumption that they define inverse 
operations. 

One can immediately notice that our free algebras provide suf- 
ficient conditions under which this assumption is enforced. This 
suggests the possibility that such generalized constructor/decon- 
structor pairs could provide the combined benefits of pattern match- 
ing and data abstraction, with the implication that direct syntactic 
support for such constructs can bring significant expressiveness to 
functional programming languages. 

3.1 Generalized Constructors with apply /unapply in Scala 

Besides supporting case classes and case objects that are 
used (among other things) to implement pattern matching, Scala's 
apply and unapply methods [5, 10] allow definition of cus- 
tomized constructors and destructors (called extractors in Scala). 

We will next describe how arithmetic operations with our AlgT 
terms, represented as ordered rooted binary trees, can benefit from 
the use such "generalized constructors". 

Our AlgT free algebra will correspond in Scala to a case 
object / case class definition, combined with a mechanism 
to share actual code, encapsulated in the AlgT trait. 

case object T extends AlgT 

case class C(l: AlgT, r: AlgT) extends AlgT 

trait AlgT { 

def s(z: AlgT): AlgT = z match { 
case T C(T, T) 

case C(T, y) => d(s(y)) 
case z => C(T, h(z)) 

} 

def p(z: AlgT) : AlgT = z match { 



case C(T, 


T) 


=> T 


case C(T, 


y) 


d(y) 


case z 




=> C(T, p(h(z))) 


} 







Note the predecessor function called p and our auxiliary functions 
named d (which "doubles" its input, assumed different from T) and 
h (which "halves" its input, assumed "even" and different from T). 



def d(z: AlgT) : AlgT = z match { 
case C(x, y) => C(s(x), y) 

} 

def h(z: AlgT) : AlgT = z match { 
case C(x, y) => C(p(x), y) 

} 

} 

We will define our generalized constructor/destructor S rep- 
resenting the successor function and predecessor function on 
rooted ordered binary trees of type AlgT by providing apply and 
unapply methods expressed in terms of our "real" constructors T 
and C and the actual algorithms defined in the (shared) trait AlgT. 

object S extends AlgT { 
def apply (x: AlgT) = s(x) 

def unapply (x: AlgT) = x match { 
case C(_, _) =>■ Some(p(x)) 
case T => None 

} 

} 

The definition of the generalized constructor/destructor D repre- 
senting double / half is similar. Note the use of the method d defined 
in the trait AlgT. 

object D extends AlgT { 
def apply (x: AlgT) = d(x) 

def unapply (x: AlgT) = x match { 
case C(C(_, _) , _) => Some(h(x)) 
case _ None 

} 

} 

The definition of the generalized constructor/destructor D can be 
seen as corresponding to \x.2x + 1 and its inverse. 

object extends AlgT { 

def apply(x: AlgT) = C(T, x) 

def unapply (x: AlgT) = x match { 
case C(T, b) => Some(b) 
case _ => None 

} 

} 

The definition of the generalized constructor/destructor D can be 
seen as corresponding to \x.2x + 2 and its inverse. Note the use 
of the generalized constructors S, D and 0, both on the left and 
right side of match statements, illustrating their usefulness both 
as constructors and as extractors. 

object I extends AlgT { 

def apply(x: AlgT) = S(0(x)) 

def unapply (x: AlgT) = x match { 
case D(a) => Some (p (a)) 
case _ => None 

} 

} 



3.2 A Scala-based Natural Number Arithmetic Package 
using AlgT Terms 

We will now illustrate how the use of generalized constructors helps 
writing a fairly complete set of arithmetic algorithms on therms 
of AlgT seen as natural numbers. For comparison purposes, the 
reader might want to look at the Haskell code in [17] where simi- 
lar algorithms are expressed using a type class-based mechanism. 
However, while the use of type classes comes with the benefits of 
data abstraction it needs separate functions for constructing, de- 
constructing and recognizing terms to express the equivalent of the 
generalized constructors used here. 

We start with a comparison function returning LT , EQ , GT and 
supporting a total order relation on AlgT, isomorphic to the one 
on N. Note here the use of the generalized constructors and I 
providing a view of the terms of AlgT as terms of the free algebra 
BinT. 



trait Tcompute extends AlgT { 




def cmp(u: AlgT, v: AlgT): Int = (u, v) 


match { 


case (T, T) => EQ 




case (T, _) => LT 




case (_, T) GT 




case (0(x), 0(y)) cmp(x, y) 




case (I(x), I(y)) cmp(x, y) 




case (0(x), I(y)) =>■ strengthen(cmp(x 


, y), LT) 


case (I(x), 0(y)) strengthen(cmp(x 

} 


, y), GT) 


val LT = -1 




val EQ = 




val GT = 1 




private def strengthen (rel : Int, from: 


Int) = 


rel match { 




case EQ =>■ from 




case =$> rel 




} 





Addition is expressed compactly in terms of the generalized con- 
structors 0, I and S. 



def add(u: AlgT, v: AlgT): AlgT = (u, v) match { 
case (T, y) =>■ y 

case (x, T) => x 

case (0(x), 0(y)) => I(add(x, y)) 
case (0(x), I(y)) => 0(S(add(x, y))) 
case (I(x), 0(y)) => 0(S(add(x, y))) 
case (I(x), I(y)) => I(S(add(x, y))) 

} 

The definition of subtraction is similar, except that the code of the 
predecessor function p is conveniently inherited directly from the 
trait AlgT, given that the trait Tcompute extends it. 

def sub(u: AlgT, v: AlgT): AlgT = (u, v) match { 
case (x, T) => x 

case (0(x), 0(y)) =>■ p(0(sub(x, y))) 
case (0(x), I(y)) => p(p(0(sub(x, y)))) 
case (I(x), 0(y)) => 0(sub(x, y)) 
case (I(x), I(y)) =>■ p(0(sub(x, y))) 

} 

The multiplication operation is similar to the Haskell code in sec- 
tion ??, except for the use of the generalized constructor 0. 

def multiply(u: AlgT, v: AlgT): AlgT = (u, v) match { 
case (T, _) =^ T 
case (_, T) T 

case (C(hx, tx) , C(hy, ty)) =>■ { 
val v = add(tx, ty) 
val z = p(0(multiply(tx, ty))) 
C(add(hx, hy) , add(v, z)) 



} 

} 

Similarly, a constant time complexity definition is given here for 
the exponent of 2 operation, by using the "real" constructor C. 

def exp2(x: AlgT) = C(x, T) 

The power operation pow takes advantage of the generalized con- 
structors and I on the left side of a case statement through the 
AlgB view of AlgT. 

def pow(u: AlgT, v: AlgT): AlgT = (u, v) match { 
case (_, T) => C(T, T) 

case (x, 0(y)) =$* multiply(x, pow(multiply(x, x) , y) ) 
case (x, I(y)) => { 

val xx = multiply(x, x) 

multiply(xx, pow(xx, y) ) 

} 

} 

Efficient division with remainder is a slightly more complex algo- 
rithm, where we take advantage of generalized constructors, direct 
inheritance from trait AlgT as well as number of previously defined 
functions: 

def div_and_rem(x: AlgT, y: AlgT): (AlgT, AlgT) = 
if (cmp(x, y) = LT) (T, x) 
else if (T = y) null // division by zero 
else { 

def try_to_double(x:AlgT, y:AlgT, k:AlgT): AlgT = 
if (cmp(x, y) = LT) p(k) 
else try_to_double(x, D(y), S(k)) 

def divstep(n: AlgT, m: AlgT): (AlgT, AlgT) = { 
val q = try_to_double(n, m, T) 
val p = multiply (exp2(q) , m) 
(q, sub(n, p)) 

} 

val (qt, rm) = divstep(x, y) 
val (z, r) = div_and_rem(rm, y) 
val dv = add(exp2(qt) , z) 
(dv, r) 

} 

Division and reminder can be separated using Scala's projection 
functions: 

def divide (x: AlgT, y: AlgT) = div_and_rem(x, y)._l 
def reminder(x: AlgT, y: AlgT) = div_and_rem(x, y) . _2 

Finally, the greatest common divisor gcd and the least common 
multiplier 1cm are defined as follows: 

def gcd(x: AlgT, y: AlgT): AlgT = 

if (y = T) x else gcd(y, reminder(x, y)) 

def lcm(x: AlgT, y: AlgT): AlgT = 
multiply(divide(x, gcd(x, y)), y) 

} 

The trait Tconvert implements efficiently conversion to/from 
Scala's Biglnt arbitrary size integers using bit-level operations 
corresponding to power of 2 and recognition of odd and even natu- 
ral numbers. The function f romN builds an AlgT tree representation 
equivalent to a Biglnt. 



trait Tconvert { 




def fromN(i: Biglnt): AlgT 


= { 


def oddN(i: Biglnt) = 




i.testBit(O) 




def evenN(i: Biglnt) = 




i != Biglnt (0) && !i. 


testBit(O) 



def hN(x: Biglnt) : Biglnt = 
if (oddN(x)) 

Biglnt (0) 
else 

Biglnt(l) + hN(x » 1) 

def tN(x: Biglnt): Biglnt = 
if (oddN(x)) 

(x - Biglnt (1)) » 1 
else 

tN(x » 1) 



if (0 = i) T 

else C(fromN(hN(i)) , 



fromN(tN(i))) 



} 



The function toN converts an AlgT tree representation to a Biglnt. 

def toN(z: AlgT): Biglnt = z match { 
case T =>■ 
case C(x, y) => 

(Biglnt(l) « toN(x) .intValueO) * 
(Biglnt (2) * toN(y) + 1) 

} 

} 

Note that for both these conversions we have used, for efficiency 
reasons, the "real constructors" T and C, although much simpler 
(and slower) converters can be built using either the AlgB or AlgU 
view of AlgT terms. 

The use of Scala's generalized constructors inspired by our free 
algebra isomorphisms has shown the combined flexibility of in- 
heritance as a mechanism for data abstraction and convenient pat- 
tern matching allowing the design of our algorithms in a functional 
style. The implicit use of apply and unapply methods in combina- 
tion with our simple free algebra semantics has facilitated the safe 
use of fairly complex (mutually) recursive functions in the defini- 
tion of the generalized constructors. The use of Scala's traits has 
facilitated flexible inheritance mechanisms supporting shared defi- 
nitions without any additional syntactic clutter. 

4. An Application: Rational Arithmetic in Scala 
with Calkin- Wilf Trees 

We will extend our Scala code snippet described in subsection 
3.1 to a realistic arbitrary size arithmetic package. It is somewhat 
unconventional, as it is based on the Calkin- Wilf bijection [3, 6] 
between N and the set of positive rational numbers Q + , rather than 
more typical representations like the arrays of long words used in 
Java's BigDecimal package, also adopted through a wrapper class 
with the same name by Scala (which runs on top of the Java Virtual 
Machine). 

Among its advantages, division (with non-zero) always returns a 
finitely represented rational and "no bit is lost" in the representation 
as canonical rational numbers with co-prime numerator/denomina- 
tor pairs are bijectively mapped to natural numbers. Our approach 
emphasizes the fact that a mathematical concept defined tradition- 
ally through equivalence classes and quotients, can be expressed 
entirely in terms of a free algebra-based mechanism. 

The trait Q representing our rational number data type contains 
distinct constructors for positive (P), negative numbers (M) and zero 
(Z). 

trait Q extends Qcode 

case object Z extends Q 

case class P(x: (AlgT, AlgT)) extends Q 

case class M(x: (AlgT, AlgT)) extends Q 



The actual code will be shared through the trait Qcode that also 
mixes-in functionality from the natural number operations defined 
in the traits Tcompute and Tconvert. 

We start with a type definition for ordered pairs of natural num- 
bers PQ represented as terms of AlgT and the conversion function 
to a conventional fraction represented as an ordered pair of Biglnt 
objects. The conversion function toFraq uses the AlgT to Biglnt 
converter toN. 

trait Qcode extends Tcompute with Tconvert { 
type PQ = (AlgT, AlgT) 

def toFraqO : (Biglnt, Biglnt) = this match { 
case Z => (0, 1) 

case M((a, b)) => (-(toN(a)), toN(b)) 
case P((a, b)) => (toN(a), toN(b)) 

} 

The function t2pq splits its argument u seen as a natural number 
into its corresponding Calkin- Wilf rational, represented as a pair of 
positive natural numbers of type PQ. Note the use of our generalized 
constructors and I distinguishing between odd and even numbers. 
The algorithm uses an encoding of the path in the Calkin- Wilf tree 
as a member of AlgB, where is interpreted as a command to take 
the left branch and I is interpreted as a command to take the right 
branch at a node of the Calkin- Wilf tree (shown in Fig. 1, for a few 
small positive rationals, represented as conventional fractions). 

I/I 

/ \ 



1/2 



2/1 



3/1 



1/3 3/2 2/3 

/. i\ t\ J\ 

1/4 4/3 3/5 5/2 2/5 5/3 3/4 4/1 
Figure 1. The Calkin- Wilf Tree 



def t2pq(u: 
case T => 
case 0(n) => 
val (x, y) 
(x, add(x, 

} 

case I(n) => 
val (x 
(add(x 

} 

} 



AlgT) : 
(S(T), 
=S> { 



PQ = u match { 

som 



= t2pq(n) 

y)) 



y) 
y). 



= t2pq(n) 

y) 



The function pq2t fuses back into a "natural number" represented 
as a term of AlgT, corresponding to the path in the Calkin- Wilf tree, 
a pair of co-prime natural numbers representing the (numerator, 
denominator) pair defining a positive rational number. 



def pq2t(uv: PQ) : AlgT = uv match 


{ 


case (0(T) , 0(T)) T 




case (a, b) =>■ 




cmp(a, b) match { 




case GT => I(pq2t(sub(a, b) . 


, b)) 


case LT 0(pq2t(a, sub(b, 


a))) 


} 




} 





This brings us to the definition of the bijection between signed ra- 
tionals and terms seen through the use of our generalized construc- 
tors and I as terms of AlgB representing natural numbers. 



def fromTCt 


1 AlgT) : Q = t 


match { 




case T 


=> Z // zero - 


-> zero 




case D(x) 


=> M(t2pq(x)) 


// odd -! 


■ negative 


case I(x) 


P(t2pq(x)) 


/ / even - 


-> positive 


} 



Its inverse from signed rationals to terms of AlgT, seen as natural 
numbers, proceeds by case analysis on the Q data type. Note that 
positive sign is encoded by mapping to even naturals and negative 
sign is encoded by mapping to odd naturals. 



def toT(q: q) : 


AlgT = q match { 




case Z => 


T // zero — > zero 




case M(x) =>■ 


0(pq2t(x)) // negative - 


-> odd 


case P(x) =>■ 

} 


I(pq2t(x)) // positive - 


-> even 



The bijection between Scala's Biglnt, seen as a natural num- 
ber type and signed rationals, is defined as the pair of functions 
rat2nat and nat2rat 



def nat2rat(n: Biglnt): Q = f romT(f romN(n) ) 
def rat2nat(q: Q) : Biglnt = toN(toT(q) ) 



Next we define a simplifier of positive fractions represented as a 
pair, to facilitate arithmetic operations on our rationals. 



def pqsimpl(xy: 


PQ) = { 


val x = xy._l 




val y = xy._2 




val z = gcd(x, 


y) 


(divide(x, z) , 


divide(y, z)) 


> 





We also use our simplifier to import and export non-canonically 
represented rationals represented as Biglnt pairs. 



def fraq2pq(nd: (Biglnt, Biglnt)): Pq = 
pqsimpl((fromN(nd._l) , f romN(nd. _2) ) ) 

def pq2fraq(nd: Pq) : (Biglnt, Biglnt) = 
(toN(nd._l), toN(nd._2)) 

We are now ready for our arithmetic operations. The template func- 
tion pqop, parameterized by a function f , will be shared between 
addition and subtraction. Note that it also involves simplification, 
to ensure that the results are in a canonical co-prime numerator/de- 
nominator form. 

def pqop(f: (AlgT, AlgT) => AlgT, xy:Pq, uv:Pq):Pq = { 
val (x, y) = xy 
val (u, v) = uv 
val z = gcd(y, v) 
val yl = divide (y, z) 
val vl = divide (v, z) 

val num = f (multiply (x , vl) , multiply(u, yl)) 
val den = multiply (z, multiply(yl, vl)) 
pqsimpl ( (num, den)) 

} 

We can use it to define addition and subtraction of positive rationals 
by simply instantiating our function parameter f to add and sub 
operating on terms of AlgT. 

def pqadd(a: Pq, b: Pq) = pqop(add, a, b) 

def pqsub(a: Pq, b: Pq) = pqop(sub, a, b) 

The comparison operation providing a total ordering of Q + relies 
on the function cmp comparing terms of AlgT seen as natural 
numbers. 

def pqcmp(xy: Pq , uv: Pq) = { 
val (x, y) = xy 



val (u, v) = uv 

cmp (multiply (x, v) , multiply(y, u) ) 

} 

Multiplication, inverse and division on Q + are defined as usual. 

def pqmultiply(a: Pq, b: Pq) = 

pqs impl (multiply ( a. _1 , b._l), multiply (a. _2 , b._2)) 

def pqinverse(a: Pq) = (a._2, a._l) 

def pqdivide(a: Pq, b: Pq) = pqmultiply (a, pqinverse (b) ) 

We are ready to define arithmetic operations on the set of signed 
rationals Q, by case analysis on their sign. We start with the oppo- 
site of a rational. 

def ropposite(x: q) = x match { 
case Z Z 
case M(a) => P(a) 
case P(a) =4> M(a) 

} 

Addition is defined by case analysis on the sign and calls to the 
appropriate operations on positive rationals. 

def radd(a: q, b: q): q = (a, b) match { 
case (Z, y) => y 

case (M(x), M(y)) M(pqadd(x, y) ) 
case (P(x), P(y)) P(pqadd(x, y)) 
case (P(x), M(y)) => pqcmp(x, y) match { 

case LT =>■ M(pqsub(y, x)) 

case Eq Z 

case GT =>■ P(pqsub(x, y)) 

} 

case (M(x), P(y)) =>■ ropposite(radd(P(x) , M(y))) 

> 

Subtraction is defined similarly. 

def rsub(a: q, b: q) = radd(a, ropposite(b) ) 

def rmultiply(a: q, b: q): q = (a, b) match { 

case (Z, _) => Z 

case (_, Z) => Z 

case (M(x) , M(y)) P(pqmultiply(x, y) ) 

case (M(x), P(y)) => M(pqmultiply(x, y) ) 

case (P(x), M(y)) => M(pqmultiply(x, y) ) 

case (P(x), P(y)) => P(pqmultiply(x, y) ) 

} 

Finally we define the inverse on non-zero rationals 

def rinverse(a: Q) = a match { 
case M(x) =>■ M (pqinverse (x) ) 
case P(x) => P (pqinverse (x) ) 

} 

and use it to derive from it a division operation on Q 

def rdivide(a: q, b: q) = 
rmultiply(a, rinverse(b)) 

} 

These operations conclude the trait Qcode. While this complete 
arithmetic package was built mostly as a proof of concept for the 
expressiveness of our free algebra based approach on progressively 
more interesting mathematical objects, future work is planned for 
turning this package into a practical tool. A first observation toward 
this end is that, like in the case of Java's Biglntegers or the C- 
based GMP package, one needs to use a hybrid approach, taking 
advantage of actual machine words (64 bits at this point), to store 
and operate on numbers that fit in a machine word. 



5. Related Work 

Numeration systems on regular languages have been studied re- 
cently, e.g. in [12] and specific instances of them are also known as 
bijective base-k numbers [20]. Arithmetic packages similar to AlgU 
and AlgB are part of libraries of proof assistants like Coq [9] and 
the corresponding regular languages have been used as a basis of 
decidable arithmetic systems like (W)S1S [2] and (W)S2S [11]. 

Arithmetic computations based on more complex recursive data 
types like the free magma of binary trees (essentially isomorphic 
to the context-free language of balanced parentheses) are described 
in [17] and [16], where they are seen as Godel's System T types, 
as well as combinator application trees. In [15] a type class mecha- 
nism is used to express computations on hereditarily finite sets and 
hereditarily finite functions. However, none of these papers pro- 
vides proofs of the properties of the underlying free algebras or 
uses mechanisms similar to the generalized constructors described 
in this paper. 

A very nice functional pearl [6] has explored in the past (us- 
ing Haskell code) algorithms related to the Calkin- Wilf bijection 
[3]. While using the same underlying mathematics, our Scala-based 
package works on terms of the AlgT free algebra rather than con- 
ventional numbers, and provides a complete package of arbitrary 
size rational arithmetic operations taking advantage of our general- 
ized constructors. 

6. Conclusion 

We have shown that free algebras corresponding to some basic 
data types in programming languages can be used for arithmetic 
computations isomorphic to the usual operations on N. 

As a new theoretical contribution, we have worked-out details of 
proofs, based only on elementary mathematics, of essential proper- 
ties of the mutually recursive successor and predecessor functions, 
on the free algebra AlgT of ordered rooted binary trees. 

A concept of generalized constructor, for which we have found 
simple implementations in Scala, has been introduced. By work- 
ing in synergy with our free algebra isomorphisms we have de- 
scribed, using language constructs like Scala's apply / unapply, 
simple and safe means to combine data abstraction and pattern 
matching in modern-day functional and object oriented languages. 

As a new practical contribution, a complete arbitrary size signed 
rational number package written in Scala has been derived working 
with terms of the AlgT free algebra of rooted ordered binary trees 
with empty leaves. 

Future work is planned to investigate possible practical applica- 
tions of our algorithms to symbolic and/or arbitrary length integer 
arithmetic packages and to parallel execution of arithmetic compu- 
tations on AlgT. 

The code snippet showing the use of Scala's apply and unapply 
methods to support generalized constructors as well as the arith- 
metic on rationals is available as a separate file at http://######. 
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